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ABSTRACT 

We use QCD sum-rules to study the decays — > 777 and — > 7/7, obtaining B((f> — > 
777) = (1.15±0.2) 10~ 2 and 13 ((f) — ► 7/7) = (1.18±0.4) 10~ 4 , in very good agreement 
with existing experimental data. We also discuss the issue of 77—77' mixing, predicting 
the 77 and 77' decay constants in a mixing scheme in the quark-flavour basis. 



1 Introduction 



The reasons for interest in (f> radiative decays are manifold. These decays yield 
information on low energy hadron physics, providing insight into the controversial 
issue of 77 — 77' mixing. From the experimental side, new data are expected from the 
KLOE experiment at the DA$NE 0-factory 0) , where a large sample of (f> decays will 
be collected, which will considerably improve the presently available statistics @) . 

In the following we study — > 777, — > r] >r y decays using QCD sum-rules ® . 
After surveying 77 — 77' mixing, we consider the coupling of the strange pseudoscalar 
current to 77 and 77', preliminary to our analysis of decays. Though we work with 
interpolating currents defined in the quark flavour basis, our results do not depend 
on any mixing scheme, providing mixing scheme independent QCD predictions. 
*Presented by F. De Fazio 



We also consider the coupling of 77 and 77' to the strange and non-strange axial 
currents, identifying the results with the decay constants in the flavour basis mixing 
scheme and estimating the mixing parameters. Finally, we give conclusions. 

2 On 7] — i]' Mixing 

77 — 77' mixing is a much debated subject. The once conventional description intro- 
duced a single mixing angle in the octet-singlet flavour basis S). More recently, it 
has been shown l§) that a proper treatment requires two angles with a redefinition 
of the particle decay constants. An equivalent description adopts the quark-flavour 
basis instead of the octet-singlet one ^ , in which such constants are defined as: 

<0|J£jP(p)>=i/&v (a = q,s; P = 77,77'), (1) 

with J| M = -^(U'j^u + dj^d), Jf^ = 57^753. The decay constants are written 
according to the following mixing pattern: 

f% = f q cos 0<? fr, = ~fs sin0 s 

f} = fq sin 0? f%> = fs cos 0s • ( 2 ) 

Feldmann HP has shown that, since \<f> 8 — (f) q \/((f> s + <f> q ) <C 1, this mixing framework 
is simpler, being specified quite accurately in terms of a single angle = (p q = <p s . 

3 77 and 7/ couplings to the pseudoscalar current 

Let us define: < 0|s275s|77 >= A, and compute this quantity using QCD sum-rules 
starting from the two-point correlator: 

T A {q 2 ) = i J d A xe iq - x < 0\T[J°{x) J 5 st (0)]|0 > (3) 

where J| = si^^s. The correlator @ is given by the dispersive representation: 

1 f°° o{ s) 

T A (q 2 ) = - ds rv \ + subtractions . (4) 

7T Hml S — q 2 

For low values of s, p(s) receives contribution from the coupling of the 77 to the 
pseudoscalar current. Hence we can write (we discuss later possible subtractions): 

A 2 1 r°° n hadf„\ 

T A {q 2 ) = ^- 2 + - ds P -^4, (5) 
m z — q z 7r J s s — q z 



assuming that higher contributions start from an effective threshold Sq. T^(g 2 ) 
can also be computed in QCD by expanding the T-product in ([|) by an Operator 
Product Expansion as the sum of a perturbative term plus non-perturbative ones 
proportional to vacuum condensates. Dispersively writing the perturbative term 
too, we have: 

1 roc rft CD (H) 

T% CD {q 2 ) = - ds £ V + d 3 <ss> +d 5 < sgaGs > +... . (6) 

7T J4ml S — q 2 

pQCD anc j ^ ^ are com puted in QCD. Now we implement quark-hadron duality, 
assuming that p had (s) and p QCD (s) give the same result when integrated above s . 
This leads to the sum-rule: 

a 2 i r° p QCD (s) 

— _ = - / ds hr + d s < ss > +d 5 < sgaGs > +... (7) 

m 7) ~ Q n J&ml s — q 2 

Now we apply to both sides of (0) a Borel transform, defined as 

B[f(Q 2 )} = Um Q2 _ >n _^ =M2 J^yi-Qr (i^JfiQ 2 ) > (8) 

where / is a generic function of Q 2 = —q 2 and M 2 is known as the Borel param- 
eter. This operation improves the convergence of the sum rule and eliminates the 
contribution of subtraction terms in (HI). The final sum- rule reads: 



A 2 e m 1 = -—r / ds s \l 1 -e ■ 
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< > 1 ~ 77? + 771 + 77? < sptrGs > 1 



• (9) 



M 2 \P) M 2 \ 2M 2 J_ 

We use < sgaGs >= 0.8 GeV 2 < ss >, < ss >= 0.8 <qq>,< qq >= (-0.24) 3 
GeV 3 , m v = 0.548 GeV. We vary the strange quark mass in the range: m s = 
0.125-0.140 GeV0) and s below the rj' pole between 0.9 2 -0.95 2 GeV 2 . Since M 2 
is an unphysical parameter, we look for a range of its values ("stability window") 
where the sum-rule is almost independent of M 2 . We fix M 2 in [0.8, 1] GeV 2 and, 
considering the uncertainty on m s and on sq, we get @): 

\A\ = (0.115 ± 0.004) GeV 2 . (10) 

Let us now consider: < 0| sry 5 s|?/ >= A'. An analogous calculation gives: 



(A') 2 e-l^ + A 2 e~^ = -_ / " dssJl- 
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where we have raised the threshold up to s' as to pick up the rj pole too. Using 
m v , = 0.958 GeV and in the stability window [1.2, 2] GeV 2 for M 2 , we obtain D: 

\A'\ = (0.151 ±0.015) GeV 2 . (12) 

Though we cannot actually establish the sign of A, A' , we assume that A ■ A' > 0. 



4 Radiative — > 777 and — > 7/7 decays 

Let us define: < 77(52) |s7^l <t>(<lh £i) >= ^V) e^ 5 (9i)a(92)/3(ei)5, (g = gi-ft)- In 
order to compute — > 777 decay, we need the coupling g = — |F(0), obtained for a 
real photon coupling to a strange quark. We consider: 

IV = z 2 | d 4 x A e -^-+^ < 0\T[J°(y)J v (0)J^x)}\0 >= II(g 2 , g 2 , q 2 ) e^q?q§ 

(13) 

with J| defined above and J v = s^ v s. The sum-rule is built up for II(g 2 , g 2 , g 2 ). 
First we write n(g 2 , q 2 , q 2 ) according to a dispersion relation in the variables qf, q\: 

ll V 2 ) = -J ds, I ds 2 ( gCfffij?) . (14) 

Now the spectral function contains, for low values of s\, S2, a double 5— function 
corresponding to the transition — > 77. Extracting this contribution, we derive the 
sum-rule (after a double Borel transform in the variables — qf, — <7 2 ): 
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The integration domain over Si,s 2 (specified in 0/) depends on q 2 . We compute 
F(q 2 ) for negative values of q 2 , and then extrapolate the result to q 2 = [J Since we 
only know the magnitude of A, it is ^(g 2 )! that is determined. We use: m^ = 1.02 
GeV, = 0.234 GeV (obtained from the experimental datum on —> e + e~ 0)). We 
use two values of the threshold: S01 = 1.8, 1.9 GeV 2 ; the 77 threshold is chosen as 
in section 3. The extrapolation to q 2 = gives \g\ = F(0)/3 = (0.66 ± 0.06) GeV -1 ; 

this way, we could perform a double Borel transform in the two variables Q\ = —q\ and 
Q\ = —q%, which allows us to remove single poles in the Si and S2 channels ("parasitic" terms). 



the uncertainty is obtained varying all the input parameters in the sum rule. From 
this result we compute T(0 — > 777) and, using T(0) = 4.43 MeV EP, find H); 

B(<j> -> 777) = (1.15 ± 0.2)% , (16) 

in agreement with the experimental result: B(<j) — > 777) = (1.18 ± 0.03 ± 0.06)% @). 
Applying the same analysis to the rj mode gives \g'\ = F'(0)/3 = (1.0 ±0.2) GeV -1 
and ID 

B((f> -> 7/7) = (1.18 ±0.4) 10" 4 , (17) 

which agrees with the experimental datum 23(0 — > 7/7) = (0.82lojg ±0.11) 10~ 4 ©. 
The results fllBf ) and (|j]) are independent of any mixing scheme for the 77 and 77'. 
Nevertheless, adopting the mixing scheme in the flavour basis, one gets: 

R = ^nl (^$)\ &n ^, (18) 

and hence cj) s = (34±g)°. The experimental ratio would give: S = (39.0±5'|)°. For 
a comprehensive survey of other results, we refer to S) . 

5 77 and 77' couplings to axial currents 

Now we investigate the 77, rf decay constants in both the strange and non-strange 
sectors in order to deduce the mixing pattern. We begin by considering the correlator 

IV(g 2 )=* |d 4 xe^<0|T[J|^)J|,(0)]|0> . (19) 

Following the procedure already outlined above, we obtain the sum-rule: 
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In the stability window [2, 3.5] GeV 2 for M 2 we get B> f° = (0.13 ± 0.01) GeV and, 
using the same technique: = (0.12 ± 0.02) GeV, /« = (0.144 ± 0.004) GeV, 
/n' = (0-125 ± 0.015) GeV. From these results we can estimate all the mixing 
parameters defined in (@), obtaining: S = (46.6°±7°), f s = (0.178±0.004) GeV and 
<p q = (41°±4°), /, = (0.19±0.015) GeV. Our results give \<p s -(t> q \l {(t>s + <P q ) ^ 0.065, 
confirming that this ratio is much less than 1 0) . 

Finally, let us consider the matrix element: < Old^J^ J 77 >= m 2 f^. The di- 
vergence of the axial current contains the anomaly: d^J^ = d^^j^s) = 2m s si r y^s+ 
f^GG, G being the gluon field strength tensor, G its dual. Therefore: 

2m s < 0|s77 5 s|77 (/) >= f^m 2 ^ - (0 



—GG 

An 



77W) . (21) 



Exploiting the previous results and fllOD, Ql^) , we derive from (f21|): 



(0 



—GG 



7]) = (0.008 ± 0.004) GeV 3 , (0 



—GG 

47T 



7]'} = (0.072 ±0.025) GeV 3 . (22) 



47T 

Both values in (|2~2"D are close to the quark model result of Novikov et al. H3). 
6 Conclusions 

We have analysed — > rjj, <fi — > rj /r y decays using QCD sum-rules. We begin by 
presenting a preliminary calculation of the couplings of the pseudoscalar current to 
T], rj' . The results derived require no assumption about 77 — 77' mixing and are in good 
agreement with available data. However, the uncertainty in the rj' case is large, and 
so the last word is left to the experimental improvement at DA$NE, for instance. 

We have also considered 77 — rf mixing, estimating the mixing parameters 
in a quark-flavour basis scheme. The existing spread of results gives us confidence 
that new experimental information will shed light on this topic too. 
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